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ABSTRACT 


Velocity and temperature fluctuation measurements were 
made over the open ocean from instruments mounted on the 
R/V ACANIA. These data were examined to determine the 
validity of present formulations and prediction techniques. 

Values of momentum flux, lee were inferred from the 
rate of dissipation of turbulent kinetic energy, €. Dissi- 
pation values were obtained from spectral estimates and 
inner scale estimates. Values of u, were examined for 
representativeness on the basis of the constant-flux assump- 
tion and by comparisons with other studies. The vertical 
variation of the dissipation rate was examined for possible 
effects of stability and wind-wave coupling. 

The momentum flux, computed from spectra, supported the 
constant-flux assumption for neutral conditions. For 
periods of instability, the stability corrections applied 
to the vertical variation of € resulted in the proper ad- 
justments toward the predicted slope. The reductions in 
momentum transfer during periods of stable stratifications 
were consistent with wind-wave coupling effects described 
by Davidson. The shape of a spectrum of temperature fluc- 


tuations was in agreement with predictions. 
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I. INTRODUCTION 


The structure of turbulence in the marine boundary 
layer has become of increasing importance in naval systems 
design. To gain a better understanding of turbulence sta- 
tistics, such as the spectral distribution of wind fluc- 
tuations, the effects of wind-wave coupling and the 
Boca daey of present prediction expressions, personnel at 
the Naval Postgraduate School, Monterey, California have 
undertaken a study to obtain the shipboard observational 
data necessary Ore accurate examination of these prob- 
lems. The group consists of personnel from several depart- 
ments; Meteorology, Mechanical Engineering, Oceanography 
and Physics, whose primary objective is to obtain observa- 
tional data and then to evaluate and improve existing pre- 
diction techniques. 

Several observational experiments have been conducted 
and were followed by data reduction and analyses periods. 
The second observational experiment is the one considered 
in this study and was conducted from 18 September 1973 to 
22 September 1973. Meteorological measurements were made 
from the oceanographic research ship, the R/V ACANIA shown 
in huguae:l, Lun the vicinity of San Nvcolas island. Pig-— 
ure 2 illustrates the "near" ocean environment and the 
ACANIA's anchorage where the majority of measurements were 


taken. 
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Figure 1. The R/V ACANIA is operated by the Department 
of Oceanography, NPS. 
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The meteorological measurements were of three differ- 
ent types. Two of these pertained to different scales of 
Measurements and the third was of the influence of the 
moving platform on mean and turbulent measurements. The 
scales considered were the turbulent scale and that which 
is associated with the mean wind and thermal stratifica- 
tion. 

The scope of this study is confined to turbulent fluc- 
tuations of velocity and temperature, although reference 
will be made to mean wind and thermal stratification mea- 
surements. The latter are important since they are re- 
quired to relate to results from other experiments or for 
extrapolating to regions or conditions where turbulence 
measurements are not made. A nearly unique problem facing 
the Navy is that present expressions for the turbulent 
processes are empirical in nature and have only been pre- 
viously validated, or at least tested, overland. 

The purpose of this investigation is to test several 
expressions in the "near" ocean environment, interpret the 
results and evaluate existing turbulence expressions. To 
accomplish this, it was decided to concentrate the study 
On specific formulations. The formulations include (1) an 
evaluation of the constant flux layer (2) the variations 
of € (the viscous dissipation of turbulent kinetic energy) 
with height and, where data permitted, (3) the variation 


of temperature fluctuations, 
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Il. OBSERVATIONAL AND THEORETICAL CONSIDERATIONS 


A. GENERAL BACKGROUND 

This study represents one of the first attempts to 
Make simultaneous, multiple-level observations over water 
from a ship. Previous measurements have been primarily 
overland or over water from near stable platforms. A 
study which closely approximates this experiment was de- 
scribed by ue et al (1973). The measurements re- 
ported in that study were made from the Scripps Floating 
Instrument Platform (FLIP) during the Barbados Oceano- 
graphic and Meteorological Experiment (BOMEX). From FLIP, 
velocity and temperature measurements were made with a 
roving probe which cycled continuously between measurement 
stations located 2, 4, 7 and 12 meters above mean sea 
level. Results in the present study will be compared to 
those obtained by Stegen, et al, as well as to other con- 
temporary studies. 

The following paragraphs will provide a review of the 
theoretical considerations needed to interpret the data. 
PLronmthese considerations) the fini ction we locity: (al) will 
be estimated from values of the dissipation of turbulent 
kinetic energy (€) and € will be obtained from spectra. 
For this reason we must examine both the turbulent kinetic 


energy balance expression and Kolmogorov's hypotheses, 
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B. TURBULENT BALANCE EXPRESSIONS 
A general expression for the turbulent kinetic energy 


balance has the form 


























Vee du. V.V 
Ce) i Sai ite at g éy 3 gi Pigs dala 
ae | 2 ) ij ox ‘ Th Wale is cai ox Is fe) 2 )] 
V.V 
) - i. 3 
= 3x, 3 =] (1) 


If all eddy transports except those along the vertical are 
neglected, the turbulent kinetic energy balance can be ex- 


pressed as 











|= 
H- 
= 
|'d 





(2) 

The time variation of turbulent kinetic energy is thus seen 
to depend on four processes represented by the four terms 
on the right hand side of Equation (2). The first repre- 
sents production of turbulent kinetic energy by the inter- 
action of the Reynolds stress with the mean shear flow. The 
second represents turbulent kinetic energy released or 
gained due to the work of vertical motions against the 
effects of buoyancy. The third represents viscous dissipa- 
tion of turbulent kinetic energy. The final term represents 
both pressure work and the flux divergence of turbulent 


kinetic energy due to transfer by turbulent motion. 
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To obtain those expressions normally used for estimat- 





ing the shear stress or eddy flux (ae = ao) on the 


basis of turbulent kinetic energy balance, the following 
assumptions and definitions are made: 
i) assume that the turbulent kinetic energy is in 


a statistically steady state 
: 2 








Vv 
c) i a 
oe Sa = 
ii) conditions are horizontally homogeneous for the 


irregularity fields (i.e., spatial variations 
of mean fluctuations are taken to be negligible 
over distances comparable to the largest scales 
of interest in the analysis). 


P) P) 
ao (  ) Sage ( ) Ss © 
dx dx, 


iii) conditions are near neutral 


g ~ 
ew AS = 
2 3° () 
fe} 
iv) define the coordinate system so 
us, = us = 0 


The above assumptions and definitions lead to the follow- 


ing turbulent energy balance expression 








-V_V =e+—~— [V (3) 


Furthermore, if the divergence terms are negligible, 


Equation (3) becomes 








-V_V = 6 (4) 


a ae > om 
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which states that mechanical production of turbulent 
kinetic energy equals dissipation of turbulent kinetic 
energy. 


Equation (4) can be rewritten as 





2 du 
u ma—- = €E (5) 
* 
dx, 
2 yc : 
where u, = “Viv, é For neutral conditions in a constant- 


stress boundary layer the following is an accepted expres- 


sion for the wind shear 


Gy 


a (6) 


where k is Von Karman's constant and z is the height. 
Combining Equations (5) and (6) yields the following esti- 


mate for u, based on measurements of € at known height z. 


Bin = heaiey 27* 


(7) 


Equation (7) can be used to both estimate the value of the 
momentum flux, ued and to examine if it is constant with 
respect to height, the constant-flux assumption. 

Table I provides the momentum flux calculations pre- 
sented by Stegen, et al with € values derived from time 
derivatives of the velocity fluctuations. These values 
will be compared to those obtained from shipboard measure- 


ments to establish the validity of the expressions and 


assumptions used in this study. 
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TABLE I 


Results of Momentum-Flux @CAlemiarione.” 


Height kze u, u, (corrected) 
(m) (cm3sec73) (cm sec” 1) (cm sec-l) 
11.69 4.02x10> 15.9 14.0 

6.70 BOGRNO® 19.0 dos 

3.89 Gesiox lon 18.5 17.4 

2.39 6.03x102 18.2 17.4 


LE ReOSMILER taken from Stegen, et al (1973) 
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Equation (7) can be rewritten as 
3 
Ue E Ribs Ase Tey eh, /k (8) 


Therefore a plot of Ln z versus Ln € should yield a 
straight line with -l1 slope. The intercept of this line 
with the Ln z axis, z = 1, will yield the estimate of 


u A plot of this shown by Stegen, et al appears in 


*x° 
Figure 3. Deviations from this particular relation could 
be interpreted for possible deviation from the constant- 


flux assumption or for influence of wind-wave coupling 


or non-neutral stratifications. 


Cc. SIMILARITY EXPRESSIONS 

The parameter € can be estimated on the basis of 
similarity expressions for the inertial subrange. Accord- 
ing to Kolmogorov's second hypothesis, the inertial sub- 
range is a region (in wave number space) where there is 
negligible dissipation of energy and where inertial energy 
transfer to higher wave numbers is the dominate feature. 

Figure 4 from Lumley and Panofsky (1964) depicts the 
inertial subrange considering spectral energy transfers. 

Based on a dimensional analysis Kolmogorov's second 


hypothesis yields the following functional relation 


HU ences ae 


where E(k) represents the three-dimensional energy density, 


k is the wave number and c iS an empirical constant. 
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Figure 3. Measured Values of € vs z, the Height Above 
the Mean Sea Surface (after Stegen, et al). 
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Figure 4. A Schematic Drawing Showing Spectral Transfer 
-and Illustrating (A) Energy Containing Region 
(B) Inertial Subrange (C) Dissipation Range 
and (D) Universal Equilibrium Range. 
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Unfortunately, E(k) is impossible to measure. It can be 
shown, however, that the one-dimensional energy density, 


$4, (kK), which can be measured, also has the form 


2735 57/8 


(kK) =a (9) 


ia 


where Q@ = constant. 

Furthermore, measurements at the scales of interest 
have to be made at a fixed point in the flow and the spec- 
tra therefore are obtained for temporal frequency, f. In 
order to use Equation (9), time and space scales must be 
related using Taylor's (1938) "frozen turbulence" hypo- 
thesis 


k = 20£/U 


where U is the mean wind speed, measured by a cup anemo- 
meter. The term "frozen turbulence’ implies that the 
turbulence remains unchanged during the time required for 
it to sweep past the probe. 

The following form of Equation (9) is often utilized 


and will be here, 


fo>..(£) = kob_._(k) = eieed Sea ee 


ll 11 ree 


The parameter € will be obtained from the spectra using 
Equation (10). The € values obtained from simultaneous, 
multiple level measurements will then be used to estimate 


2 ‘ 
the value of the momentum flux, u, , from Equation (7). 
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There is another way of estimating € based on the 
scale associated with the upper limit of the inertial sub- 
range. If the spectrum is in statistical equilibrium it 
is independent of the source of the turbulence and depends 
only on €, y and k. A dimensional analysis yields the 


following expression for the inner scale, be 


1 Serene 


a (11) 


where eas 271 /u) is the dissipation wave number and y 

is the kinematic viscosity of air. Pu aa eos Reus wave 
number can be obtained from spectra which show scale 
separation between pure inertial transfer (-5/3 regions) 
and the viscous dissipation region. Dissipation of turbu- 
lent kinetic energy can then be obtained from Equation (11) 


and compared to those obtained from Equation (10) to estab- 


lish greater confidence in the spectral estimate. 


D. FURTHER THEORETICAL CONSIDERATIONS 

Although the expressions previously described are based 
on accepted and consistent theories, there are several 
serious assumptions which must be considered in applying 
them to observational data. In many instances the theory 
is empirical and the formulating analyses were based on 
overland data. A recent effort in turbulence theory has 
been to prove or disprove the validity of overland formu- 
lations for non-neutral conditions or over water regimes 
where the wave influence may be important. The latter has 


been described by Davidson (1974). 
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In developing the turbulent kinetic emergy balance, 
given by Equation (3), the existence of meutral conditions 
was a critical assumption. For non-neutral stratifica- 
tions, Equation (7) would include a buoyamcy term and be 


of the form 








BS Zz z 
CP ies os Ms) = =! (12) 
or 
re 3 
ary kz et) 


where L = - Tu, /gkV,6 = Monin-Obuhkov Length; v3 = heat 
flux, 

Ri = Richardson Number, 

%, = empirical function of = ’ 

%, = empirical function of Ri. 


For the non-neutral case, the expression corresponding to 


Equation (8) becomes 


eae) 
u, 





Ln € = -Ln z+ Ln >. (Ri) + Ln (13) 


Possible deviations from the constant flux assumption or 
wind wave coupling influence exhibited by spectral results 
could be examined by testing overland expressions for $., (Ri). 
A stability correction could then be applied for periods of 
non-neutral stratifications. A possible effect of the 
instability on the predicted plots of Ln € versus Ln z is 


illustrated in Figure 5, 
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Figure 5. An Illustration of a Possible Effect of 


Stability on the Variation of € with Height. 









(Qe¥ nd 


2ep2 hoyioeewen 


ff 
Fa 
eoilserios S F 


t.e278 eldiasod «a tea nakgarseulist nh |e sor? 
| ' 16 AOLISBLIAV GHD Bu <tilbdore 4 


‘ oe 


ie —— a 


III. OBSERVATIONS 


A. THE EXPERIMENT 

The data considered in this study were obtained from 
Measurements made aboard the R/V ACANIA, which is operated 
by the Department of Oceanography of the United States 
Naval Postgraduate School. As previously stated, San 
Nicolas Island was the experimental site and the majority 
of data were collected while she ship fia at anchorage near 
the northwest tip of the island. Additional data were col- 
lected while the ship was enroute between Monterey and San 
Nicolas Island and also while the ship was transiting to 
the west side of the island. Since the data collected from 
anchorage comprise the largest percentage of successful 
runs and represents all stability conditions, data runs 
selected for analysis in this study are primarily anchorage 
data. Table II is a listing of the various data runs ex- 
amined. Figure 6 shows the ship's track in the vicinity of 
San Nicolas Island and Figure 7 shows the route to and from 


the island. 


B. DESCRIPTION OF EQUIPMENT AND MEASUREMENTS 
Figure 8 shows the sensor placement at various levels 
along with the variables measured at each level. 
Ll. Velocity Measurements 
Velocity fluctuation measurements were made at three 


levels on the forward mast and one level on a Single bow 
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TABLE II 


Data Periods Considered in This Study 


Date Time 


20 Sept. 73 1103-1138 


20 Sept. 73 2008-2018 


20 Sept. 73 2023-2053 


Digital Channel 


B WN 


NO 


Own f& W 


&m W DN FF 


Parameter 


upper 


level 


lateral accelera- 


tion 
upper 
tical 


upper 
tical 


level ver- 
acceleration 


level ver- 
acceleration 


(integrated) 





tthe subscript refers to the level at which the 
Level 1 is the uppermost 


measurement was made, 
level. 


26 





19deme1eT sense taateia eat ; a 
Pie i stitegori” 26 
- a 
2” s 
iY ; (‘waa 
, . . 
ty : atot-eoos ef) 
4 a 
«” ¢ 
' 
a” ‘ = 
’ 
»? , 
Po é 
fiwepl veyqya 8 
-Baeleone teraza! 
0e23 
“zou isvel toqan t 
rocco omen ligoon Lleol 
~teo. devel 10496 Ly 
nolitetzsieuos fond) 
(ha Jaatpedar) 
ir t CPOS~Llip OF. «i 
Lu c : 
F 
a t 
> 
P : 
io Aol te buvel ofy at #28208 Baleoadve var! 
' tejq/ @o) el | lower wie WE eden 


es 


‘ybua2 wld? wi bewabiened. sooivas : : 
ie 













eee 3 





elavat 





TABLE II (continued) 


a 


Date Time Digital Channel Parameter 


20 Sepe, TPB ASAIoAM u 
i] 
9) 
0 
3 
upper level fore-aft 
acceleration 


bm W DH F 


5 upper level lateral 
acceleration 
i] 
6 Th 
7 _ upper level vertical 
acceleration 
(integrated) 


8 upper level vertical 


acceleration 


DileSiopite a Sa OS 0.6152 1 upper level fore-aft 
acceleration 


own & WN 
4 


upper level vertical 
acceleration 


7 upper level lateral 
acceleration 


8 upper level lateral 
acceleration 
(integrated) 


ee ee ea ee eae rT 
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Figure 7. The San Nicolas Island Region. 
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mast as shown in Figure 8. Thermo-Systems (TSI) Model 
1210 probes with sliding support shields were selected 
for these measurements. The sensor consisted of a plati- 
num coated tungsten wire 0.00015 inches in diameter and 
0.060 inches long, small enough to resolve the viscous 
dissipation scale without needing to apply wire length 
corrections. The wire was fitted with plated ends for 
isolating the sensing area and thereby minimizing flow 
disturbance. 

Electronics associated with Bee probe were a TSI 
Model 1054B Linearized Anemometer, a TSI Model 1056 Variable 
Decade Module and a Model 1057 Signal Conditioner. Schema- 
tics of the anemometer circuits appear in Appendix A. The 
anemometer had a linear frequency response from DC to 10 KHz 
and the variable decade module operated with a 0-60 ohm 
range. The signal conditioner had a linear frequency 
response from DC to 400 KHz and permitted voltage suppres- 
sion in one volt steps from 0-29 volts. The voltage output 
from the conditioner was recorded on magnetic tape. 

Sensor placement required 200 foot cables instead 
of the usual 15 foot cables. This resulted in a decrease 
in system frequency response but the decrease had little 
effect in the frequency band of interest. Ten probes were 
available for the experiment and they were placed into 
four groups. Each group of probes was calibrated with 
separate cables and electronics which permitted a complete 


system calibration for each of the four sensing levels. 
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The calibration was completed with a manometer and a TSI 
Model 1125 Calibrator shown in Figures 9 and 10. Empiri- 
cal relations used in interpreting the hot wire signals 
are presented in Appendix B along with the computer code 
used to evaluate these measurements and resultant calibra- 
tion curves for each probe. 

The performance of the hot-wire system was excel- 
lent throughout the experiment. Only one probe failure 
was experienced despite continued probe exposure during 
the data gathering period. 

2. Temperature Measurements 

Temperature fluctuation measurements were made at 
two levels on the forward mast of R/V ACANIA as shown in 
Figure 8. Electronics used were TSI Model 1044 Tempera- 
ture Modules. The unit power supply consisted of six volt 
batteries. It can be utilized to make both temperature 
and AT measurements. The temperature module is designed 
for measurements of temperature variations down to 0.001°C 
at frequencies up to 1 KHz. The schematics for the Model 
1044 Temperature Module are shown in Appendix C. 

Probe failure occurred at the beginning of the 
first data run and to permit temperature measurements to 
continue the probe unit was fitted with a 60 ohm platinum 
wire sensor. 

ae Parallel Measurements 
Mean measurements of the wind speed, temperature 


and relative humidity have been described by Cavanaugh 
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(1974). Measurements of ship's motion have been described 


by Welsh (1974). 


Cc. SIGNAL PROCESSING AND RECORDING 
1. Analog Recording 

All data were recorded simultaneously in analog 
form with a 14 channel Sangamo Model 3562 FM Recorder. 
All recordings were made at a recording speed of 7 1/2 
inches per second. At this recording speed the recorder 
had a frequency range from DC to 25 KHz. 

A schematic diagram of the complete recording sys- 
tem is shown in Appendix D. The amplifier and suppressor 
combination in series with the sensor electronics and re- 
corder channels were necessary to improve the signal to 
noise ratio of the inputs. The hot-wire signals were high 
pass filtered at 2 Hz during the majority of the runs. 
Two data runs were made without the filter to permit low 
frequency correlation between the velocity signal and 
ship's motion. The results of this correlation are dis- 
cussed by Welsh (1974). The temperature signals were low 
pass filtered at 500 Hz. In addition to the signal re- 
cordings, a voice track was utilized to record observa- 
tions pertinent to the data runs. Constant monitoring of 
amplified and preamplified signals was performed with an 
oscilloscope to evaluate the signals and to insure that 
the signal would not exceed amplifier or recorder input 


limits. 
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2. Real Time Analysis 

Real time analysis was accomplished using a Spec- 
tral Dynamics Corporation (SD) Model SD301C Real Time 
Analyzer. This unde was used in conjunction with a Model 
SD309 Ensemble Averager, a Model 13116-2A X-Y Display 
Oscilloscope and an X-Y plotter to provide an instantaneous 
evaluation of the hot wire and eemperakince signals. The 
equipment allowed positive determination that measurements 
were sufficient to define the inertial subrange. A system 
operation diagram is provided in Appendix E along with 
several of the X-Y plots obtained from the real time 


analysis. 
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LV. DATA REDUCTION AND ANALYSIS 


A. PRELIMINARY ANALYSIS 

The analog tapes were played on a strip chart reconder 
to permit initial evaluations of the data. From this, 
data periods were chosen for subsequent reduction and 
analyses. This preliminary analysis disclosed a recording 
problem that had previously gone unnoticed. A periodic 
tape transport vibration had occurred due to the 7 1/2 ips 
gear drive and this obliterated the recorded signal. For- 
tunately, there were sufficient time intervals between 
fibReetiome to permit 20 minute data records to be obtained. 
It was from these data sections that the periods listed in 
Table II were obtained. 

The difficulties encountered in making the temperature 
fluctuation measurements, as described in Section III, 
were evident in the strip chart recordings. The tempera- 
ture data either failed to show sufficient signal strength 
or displayed a DC drift which rendered them impossible to 
digitize. An example of the two cases ifs shown in Figure 
Waibe 

The paucity of temperature data led to a focusing of 
attention on the velocity fluctuation measurements. The 
criteria for selection of subsequent data periods for in- 
clusion in the reduction and analysis phase consisted of: 


(1) signal strength (2) averaging time (3) DC drift rate 
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and (4) availability of mean profiles. Sufficient signal 
levels were required to minimize amplification during the 
analog-to-digital conversion process. Twenty minute data 
runs proved sufficient to permit ample averaging time for 
Fast Fourier Transform (FFT) calculations. DC drift was 
kept at a minimum to permit ease in digitization and to 
allow greater confidence in the data samples. Finally, 
data periods were selected to correspond to times when mean 
wind information was available. A typical signal meeting 


the above mentioned criteria is shown in Figure 12. 


B. ANALOG-TO-DIGITAL CONVERSION 

Analog-to-digital conversion was performed on a hybrid 
computer system located on the fifth floor of Spanagel 
Hall at NPS. The hybrid system consists of an analog com- 
puter, COMCOR Ci 5000, interfaced electrically with a 
digital computer, XDS 9300. A detailed discussion of the 
operation of this computer facility was given by McKendrick 
(1972) which should be referred to for additional informa- 
tion. A block diagram, providedin Figure 13, illustrates 


the analog-to-digital conversion sequence followed in this 


study. 


Cc. SEVEN-TRACK TO NINE-TRACK CONVERSION 
The analog-to-digital conversion process generated 
seven-track octal base data samples written on digital tape. 


The computer used to analyze the digital data was an IBM 
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Bigure 12. Typical Velocity Signal. 
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360/67 computer. This computer utilizes nine-track hexi- 
decimal based digital samples so a conversion between the 
seven-track and nine-track systems was necessary. A dis- 
cussion of theprocedures necessary for making the proper 
conversion to the nine-track hexidecimal based system are 
described in Naval Postgraduate School Computer Center 
Technical Note #0211-08. The two primary programs involved 
in this conversion are the TAPEOUT and CONVERT programs. 
These programs along with a description of digital tape 


operations are discussed in Appendix F. 


D. DIGITAL ANALYSES PROGRAMS 

- Fast Fourier Transform (FTOR) and spectral analysis 
(SCOR) programs were available in the IBM 360/67 subroutine 
library. The following description of the FFT progran 
package was taken from McKendrick. The programs were 
written by J. F. Garett and J. R. Wilson while students 
at the Institute of Oceanography at the University of 
British Columbia. The "OCEAN" subroutines of the UBCFTOR 
program instructs the IBM 360 to read the input tape, one 
block at a time, compute the FFT of the data block and 
then store the resulting Fourier coefficients on the out- 
put tape. The computer printout for the FTOR program is 
a listing of each digital channel and the Fourier coeffi- 
cients produced. An example of this output is shown in 
Figure 14, 

With these coefficients, the UBCSCOR program is used 


to compute the spectrum for a given channel, or the cross 
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FTOR Output. 


Figure 14. 
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spectrum between a pair of channels, depending on what is 
selected by the user. The output of the SCOR program con- 
sists of a tabular computer printout, as shown in Figure 

15, and a graphical plot of the spectrum versus frequency 

as shown in Figure 16. Appendix G of this thesis contains 
the preface comments of the FTOR and SCOR programs along 
with additional information for more efficient utilization 
of the spectral analysis package. Several difficulties 

were encountered during utilization of the analyses programs 
and consideration of these difficulties is the primary 


reason for inclusion of Appendix G. 
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SCOR Spectral Statistics. 


Figure 15. 
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SCOR Spectral Plot. 


Figure 16. 
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Vv. RESULTS 


A. ONE POINT SPECTRAL METHOD 

The procedure first selected for estimating € , and 
hence u,, was to use single point values from the spectra, 
similar to that in Figure 16. A -5/3 slope was approxi- 
Mated for the spectrum and the frequency and spectral den- 
sity, for a spectral value coincident with the line, were 
taken from computer output similar to that illustrated in 
Figure 15. Equation 10 was then used to obtain € from 
values of frequency, spectral density, wimd speed and probe 
height. The computed value of € was then used in Equation 
7 to obtain a value for the friction velocity (u,). The 
computer program used for these calculations appears in 
Appendix H. 

U, values were obtained by this method for-all data 
periods and tabulated results for one minute and 10 minute 
averages are listed in Tables III-VII. This method ap- 
peared successful for spectra such as those shown in Figure 
16, which had a definite -5/3 slope for a majority of 
points. aeegee the single point spectral method is quite 
subjective when used with spectra such as the one shown in 
Figure 17. It is easily seen that selections of unrepre- 
sentative spectral estimates could result in values of € 
which would differ significantly. Because of the subjec- 


tive manner in which -5/3 slope lines could be drawn and 
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TABLE III 


Momentum Flux Calculations for 1103-1138 20 Sept 1973 
using the One Point Spectral Method 


Height Time Average € el. 
(m) (min) (emia SCE) (cm eee mah) 
14.41 1 1.078 8.727 
10 1.319 9.129 
13.20 1 1.010 8.532 
10 1.389 9.109 
8.39 2.973 9.871 
10 2.075 8.865 
4.52 2.418 7.602 
10 2.636 7.812 
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TABLE IV 


Momentum Flux Calculations for 2008-2018 20 Sept 1973 
using the One Point Spectral Method 





Height Time Average € u, 
é 2 vi _ 
(m) (min) (cm sec =) (cm sec ) 
14.41 1 9.742 17.674 
10 8.937 17.272 
8.39 1 16.211 18.291 
10 11.679 15.769 
4.52 1 20.806 15.643 
10 16.128 14.287 
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TABLE V 


Momentum Flux Calculations for 2023-2053 20 Sept 1973 
using the One Point Spectral Method 


Height Time Average € u, 
(m) (min) (ome seer) (cm 3067 >) 
14.41 1 9,304 17.505 
10 . 95516 : 17.662 
8.39 1 14.042 16.768 
10 14.519 16.956 
4.52 1 21.934 15.829 
10 22.634 15.996 
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TABLE VI 


Momentum Flux Calculations for 2323-2341 20 Sept 1973 
using the One Point Spectral Method 


Height Time Average eae: ue 
(m) (min) (ene seen) (cm Soe 2) 

14.41 1 4.069 14.925 

10 3.543 12.688 ~ 
13.20 1 3.316 12.973 

10 2.576 11.082 
8.39 1 7.965 13.859 

10 6.265 12.813 
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TABLE VII 


Momentum Flux Calculations for 0637-0652 21 Sept 1973 
using the One Point Spectral Method 


Height Time Average € Wh. 
(m) (min) toma eee) (cm Secua) 
14.41 1 1.923 10.146 
10 1.114 8.625 
8.39 1 7.830 13.991 
10 5.660 12.386 
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Figure 17. Velocity Spectrum Showing Deviation from 
-5/3 Slope. 
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a coincident spectral value selected, there was some 


doubt expressed in results from other tham ideal spectra. 


B. INTERCEPT METHOD 

To avoid the subjectiveness of estimates of € based 
on single spectral estimates, a second method was utilized. 
This method is based on Equation 10 rewritten in the fol- 


lowing form 


2/3 


Ln ko?(k) = -2/3 Ln k + Ln Ge (14) 


For this expression a plot of Ln ko¢(k) versus Ln k would 
result in a curve with a slope of one and an intercept 
2/3 


value equal to Ln ae - From this intercept, the value 


of € can be obtained from 


Exp b,3/2 

[PER 2)°/? 2 
Qa 

where b = intercept value. 


The advantage in using the intercept method versus 
the single point spectral method is that the intercept is 
determined from an average of the points along a line de- 
fined by the -2/3 slope. Figure 18 is an example of most 
curves obtained from the application of Equation 14 to 
spectral results. Difficulty in approximating the slope 
was encountered in a few cases such as that illustrated 
in Figure 19 but in most instances the forcing of a +l 


slope permitted easier detection of spectral harmonics 
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which did not fit the -5/3 criteria. Table XI does not 
indicate a one minute spectral average since the curves 
for that time average were similar to Figure 20. It is 
clear that fitting a +l slope to this scattering of 
points would be difficult. 

The intercept method was used to calculate € and Un 
for all five data periods. The results of these calcula- 


tions are listed in Tables VIII-XII. 


C. COMPARISON OF METHODS 

A comparison of Tables III-VII and Tables VIII-XII 
indicates that somewhat higher u, values were obtained 
with the intercept method. The apparent validity of the 
higher values will be examined later. 

Considering the intercept method, a majority of cases 
yAenidiedualarng cram vicrlucism fOr Longerait ime avenages. janie 
importance of resultant higher u, values will be apparent 
when comparisons are made between shipboard results and 
results from other studies. The subjective nature of 
selecting spectral estimates in the one point spectral 
method is manifested by larger deviations in u, values 
associated with periods when spectra were not "near-per- 
fect." 

The one point spectral method was compared with the 
intercept method because of interests in designing pro- 
cedures for quick estimates of u, during measurements. 


If valid values of € could be obtained from single spectral 
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TABLE VIII 


Momentum Flux Calculations for 1103-1138 20 Sept 1973 
using the Zero Intercept Method 


Height Time Average | € u, 
(m) (min) den ec) (cm seen) 

14.41 1 1.329 9.151 
4 1.378 9.330 
10 1 23)5 9.012 
20 1.372 9.248 
13.20 1 1.445 9.139 
1.379 8.998 
20 1.503 9.260 
8.39 2.635 9.600 
4 2.337 9.224 
20 2.622 9.584 
4.52 3.086 8.233 
4 3.148 8.288 
20 3.689 8.738 
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TABLE IX 


Momentum Flux Calculations for 2008-2018 20 Sept 1973 
using the Zero Intercept Method 


Height Time Average € uy, 

(m) (min) (ema See 7) (cm See") 

14.41 1 10.288 18.012 

8.285 17.135 

10 - 10.935 : 18.474 

8.39 1 13.226 16.437 

11.828 15.988 

10 13.678 16.622 

4.52 1 17.765 14.731 

14.677 13.846 

10 18.127 14.855 
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TABLE X 


Momentum Flux Calculations for 2023-2053 20 Sept 1973 
using the Zero Intercept Method 


Height Time Average € uy 
(m) (min) deme Bee) (cm sear) 
14.41 1 12.071 19.093 
11.825 18.985. 
20 11.540 18.809 
8.39 1 16.723 17.774 
16.496 17.694 
20 16.243 17.603 
4.52 1 25.023 16.541 
3 24.492 16.432 
20 24.149 16.346 
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TABLE XI 


Momentum Flux Calculations for 2323-2341 20 Sept 1973 
using the Zero Intercept Method 


Height Time Average € u, 
(m) (min) (em? eke? (cm Seer) 
14.41 3 12.925 19.533 
10 13.564 19.849 
13.20 3 13.471 19.235 
10 14.009 19.488 
8.39 3 14.418 16.917 
10 15.025 17.151 
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TABLE XII 


Momentum Flux Calculations for 0637-0652 21 Sept 1973 
using the Zero Intercept Method 


Height Time Average € u, 
(m) (min) amie Sen) (cm Secnien) 

14.41 1 1.209 8.868 
1.955 10.407 
10 1.570 9.675 
8.39 1 5.316 NAG dbS}al 
8.451 14.158 
10 7.643 13.691 
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estimates then the structure function measured may be 
valid. The latter has been shown to be a poor approach 


On the basis of these results. 


D. EVALUATION OF THE VARIATION OF U, WITH HEIGHT 

A long accepted concept in turbulence theory has been 
the existence of a constant-stress layer. This concept 
is often used to evaluate the representativeness of accumu- 
lated data (e.g., Stegen, et al) and will be applied to 
test the validity of assumptions and expressions need in 
the present study. 

Table I was compiled by Stegen, et al for conditions 
when the measured mean wind speed was 5.12 m SOOk ue This 
mean speed corresponds to the wind conditions for the 
values tabulated in Tables VII and XII. A comparison of 
results in Table I and Tables VII and XII reveals lower 
values of u, in the present study. As previously men- 
tioned, the constant, (Q), in Equation 10 is aa oo he 
derived and several values have been suggested in the 
literature. A value of @ = 0.60 was used in compiling 
the results in Table I while a value of 0.48, suggested 
by Pond, et al (1966), was used in the present study. 

UP seswlhse las ced mani haber Vili aisupporty 2 Constant— 
flux assumption as well as the validity of Equation 7 for 
u, calculations, for neutral or near neutral conditions. 


The anomalously low value of u, at 4.52 meters could be 


the result of acceleration due to the ship's bow. It is 
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seen that as the spectral averaging time increases the 
value of the flux increases. This would be consistent 
with an interpretation that increased averaging time 
would reduce the observed acceleration effects. 

Tables IX-XI correspond to periods of stable strati- 
fications while Table XII reflects a period of unstable 
stratification. It is seen from Equation 12 that neglect- 
ing diabatic effects would introduce a bias into the data. 
A stability correction could be applied to adjust the flux 


values for the more general case of stratified flow. 


E. THE VARIATION OF € WITH HEIGHT 

Figures 21-25 show the plots of Ln € versus Ln gz for 
the periods listed in Tables VIII-XII. On the basis of 
Equation 8, a plot of Ln € vs Ln z should yield a -1l 
slope for neutral conditions. Results from Stegen, et al, 
Figure 3, were in agreement with this prediction. 

In Figures 22 and 23 the slopes are less than -1l while 
in Figure 24 the slope is greater than -l. An initial 
concern when these results became known was whether possible 
attenuation of the sensors at lower levels was a factor. 
Such a concern is logical since the sensors near the sur- 
face would be likely to lose sensitivity more rapidly as 
a result of increased salt loading at lower levels. 

The data period which occurred last is the one which 
shows the greatest slope which indicates that sensor deteri- 


Oration at lower levels was not responsible for the 


65 







Sherer 2ads ad Bicow ekAt 
ghis gotpatsys beens 
sQensiage "aortsaeAzaloy7s 
=lgesdu bidate 2 =hbicaq © 
‘efftedencs $6. Goitsy & edce! 
-acelfen eeas ££ warttcupd foxy? ade 


4 
- 


aa 
1 
i 
- 
4 


segeb eos yognt S619 & esuberdn { 
wee etd 20jfhe «3 batt¢ae sd Hiuo 


pwuLs BOLZiXexrss. Fo sass is 
Tao | W 
Sak SA Siersy = oI to. atoly ? 
TOsmk Ont ao eT TA-SIEV 
i= ~~ binwhy Blue, Z 
(48 24 4) nBpas2 eo12 iy 
swusrUroes z 
eiifiw f= fats eaet =x 
Lo A _ = 
@ich wOc sed7Taiw ‘ rc 
( a f , 
i ; 
« int “ a r 
isovek t 
4 2 { 
*iaeJo youl 5 


mod. ageherodt smiz baiporwysé Inia vole yar on ade 
evebutau! ae etd to 
tara foil aro Teer irk re 
u% KOvidcedsa #dg anuhett 


S(SgGeets03 











Iin=x%r Bed 


dae? slindw 228 
32 .coltest® 
Foattse 


nao. ss 
iy « 
f 
‘ t 
Ay gaat? 
a ; aw od] 
Ss acotteLe 


Lne 


2.5 


2.0 


1.5 


1.0 


0.5 





) 0:5 1.0 5 2:0 25 


LnZ 


Figure 21. Measured Values of € vs Height for 1103-1138 
20 Sept 73. For Scaling Purposes a Factor of. 
5.0 was Subtracted from the x-Values. 
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Measured Values of € vs Height for 2008-2018 
20 Sept 73. For Scaling Purposes a Factor of 


5.0 and 1.4 were Subtracted from the X and Y 
Values, Respectively. 
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Figure 23. Measured Values of € vs Height for 2023-2053 


20 Sept 73. For Scaling Purposes a Factor of 


5.0 and 2.0 were Subtracted from the X and Y 
Values, Respectively. 
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Figure 24, Measured Values of € vs Height for 0637-0652 
21 Sept 73. For Scaling Purposes a Factor of 
6.0 was Subtracted from the X Values. 
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Figure 25. 
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Measured Values of € vs Height for 2323-2341 
20 Sept 73. For Scaling Purposes a Factor of 
5.0 and 1.4 were Subtracted from the X and Y 
Values, Respectively. 
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variation in slopes. Further support of this conclusion 
is provided by using Equation 11 as an imdependent method 
of determining the dissipation. 

For the 2023-2053 20 Sept 73 data period, the spectra 
at all levels were similar to that shown in Figure 26, 
where scale separation between the pure imertial transfer 
and the viscous dissipation Legions was discernable. The 
upper limits of the inertial subrange were determined from 
the spectra and Equation 11 was used to obtain an estimate 
of the dissipation. In Table XIII the values obtained by 
this method are compared to those obtained from spectral 
estimates. The inner scale computation is independent of 
sensor calibration since it is determined by frequency dis- 
tributions and not on the spectral value.w The apparent 
agreement between the two estimates is am in situ calibra- 
tion which indicates the loss of sensitiwity by individual 
sensors was not the cause for the deviation from the pre- 
dicted slope. 

The deviation from the predicted slope was also ex- 
amined for possible influence of non-neutral stratifica- 
tions or wind-wave coupling. Deacon (1951), from examina- 
tions of numerous vertical wind profiles, proposed the 


following expression for the wind shear 
Z 
Ne ee G ) (15) 


where £8 is a function of stability. 
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Figure 26. Twenty Minute Velocity Spectrum for 2023-2053 
20 Sept 73. 
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TABLE XIII 


A Comparison of Dissipation Obtained by Spectra (€) and 
Dissipation Obtained from the Inner Scale (e.) 

= Nhs Seo ee Ie See ee ee ee 
ee oeoa=S= = 0006090060 


Height € € ay 
(m) ee, -3 2 a3 


14.41 9.556 9.516 17.662 
8.39 14.519 14.304 16.956 
4.52 22.634 22.492 15.996 
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The variation of 8 as a function or Ri, proposed by Dea- 
con, is shown in Figure 27. In Deacon's results, the 
variation of B may be described in terms of hydrostatic 
stability. For neutral conditions, 8 = 1 so Equation 15 
then reduces to Equation 6 and the -1 slope on the Ln € 
versus Ln z plot is predicted. 

For unstable conditions 8 is greater than 1 and for 
stable conditions f§ is less than l. Deacon's results 
predict that the expected slope of a Ln € versus Ln z 
plot would deviate from the neutral case with slopes less 
than -1l for periods of stable stratifications and with 
slopes greater than -l for unstable conditions. Examina- 
tions of Figures 22-24 indicate why such considerations 
are important in interpreting the shipboard results. 
Figures 22-23 correspond to periods of stable stratifica- 
tions with the 2023-2053 20 September data period being 
the most stable. In both figures, slopes less than -l 
occur, with the more stable of the two periods showing 
the greatest deviation fein the -1 slope. The results 
are, therefore, consistent with arguments proposed by 
Deacon. Figure 24 corresponds to an unstable period and 
shows a slope greater than -l as predicted by Equation 15. 

The results were further examined with respect to pre- 
dictions from the Monin-Obuhkov similarity theory. The 
theory states that in the boundary layer there exists a 
characteristic length determined by the stability and tur- 


bulence which can be used as a scaling parameter. For 
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neutral conditions the appropriate characteristic length 
is the height above the marine boundary. 

Equations 12 and 13 show the effect of non-neutral 
considerations in developing an expression describing the 
vertical variation of €. Equation 12 indicates a relation- 
ship between the scaling parameters Ri amd L. Sucha 
relationship was described by Busch (1573) in a review of 
information available in the literature. As presented by 
Busch, the function 0, (e)= = in Equation 12 has the form 
(Qe - =) where several flux-profile relationships have 
been suggested for one Table XIV includes the formula- 
tions which Busch used to describe on along with constants 
estimated by Paulson (1970). Relationships between the 
Monin-Obuhkov length and the Richardson mumber were first 
used to obtain stability length (L) corrections for the 
shipboard results, since Richardson numbers were available 
from profile measurements described by Cawanaugh (1974). 
Values of L were then used to define (O, = = from rela- 
tionships in Table XIV. 

Figures 28 anc 29 show the effects of stability correc- 
tions for the periods with stable strattf£ications listed 
in Tables IX and X, respectively. Comparisons of the cor- 
rected and uncorrected results indicate that adjustmerts 
were toward the predicted -1l slope. The fact that the 
stability correction produces a -1l slope supports the basic 
premise of the Monin-Obuhkov theory and the validity of 


the relationships in Table XIV. 
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TABLE XIV 


Flux-Profile Relationships Describing the Stability 
Function, one 


We Stability Condition 
zZ zZ 
Bs —- > 
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Figure 28. Effect of Stability Correction for 2008-2018 
20 Sept 73. 
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Figure 29. Effect of Stability Correction for 2023-2053 
20 Sept 73. 
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The Monin-Obuhkov theory and the flux-profile relation- 
ships were also examined for the period with unstable stra- 
tification. Figure 30 shows the results of this examination. 
For -0.05 < Ri < -0.15 there is a tendency for adjustment 
toward the predicted -l slope. However, further decreases 
in the Richardson number caused an unexpected reversal of 
the adjustment. 

It is seen from Figure 30 that a value of -1.0 for the 
Richardson number results in a vertical wariation of € which 
deviates significantly from the predicted slope. The reason 
for this deviation can be explained by examining the func- 
tion, > (= - =) Figure 31, from Garratt (1972), shows a 
plot of this function using the expressions in Table XIV. 

To provide the correction necessary to obtain a -1l slope, 
(9. - =) must decrease for decreasing values of = le) aus) 
apparent from Figure 31 that this criterion is not met. On 
the basis of these results there is some evidence that the 


Gea cannot provide the adjustment 


expression g. (1-16) 
necessary tc correct the plot for conditions of large 
instability. 

The fact that the stability corrections did not com- 
pletely restore the -1l slope, even in the stable case, 
leads to examinations of the possible effects of wind-wave 
coupling. Davidson and Frank (1971) observed wave-related 


velocity fluctuations and wave-related momentum transfer 


in spectral results obtained over natural waves. Davidson 


80 


entiteiys stitorgq-xvt? add Die 


en 


7s0s644 VfoxitudO-niaok edt 


~psdn videtenu A2iw boiasd sits 


=o} 


honimane cela see 


elsandnexs divit Lo wrtlueat si3 iwote Of aiypig woods 


Oneeseuros zor PYoushusy € 2i srsiia ¢f.0= > dh > 20, Oe 
BGOBSS209 setiguut ,rsyowoli .seofe- [- EBsivotberg eda @ 
to LeexzsVox SSiJ72eeeent née S542 a0 tadmyh - $xGéus ER 


ea2 so? 7. L- u 6 aig C j nor? nese <eL 23 
eeLtw = To nattaitrae lavtria2rsy 's bo 3 rime: nets | 
eoaa0s oT -Sa0ie Gersoibuarvcd iz J aaa TH 24% 

wJawu2 O09 gir itses yt bo ne vol aida 

B Bwede ,. (ste IJ7eit3she BOT? vols _(—'< =) 8 S 
e¥aa SideT ui a ‘ ! j ity te 
seqolts I- ’ ; ivotq@ 
it H 2 seal E23 c ? 1 te wan (2 - 
a 
ao —sa8 JC i . . iq 
mit? a2 Th ; : o ttend 
PEE: eee 
1 *4eee 
‘$ilidase 
é ( i 3 i"? 
; f Tr lately 

i ’ wheal 

: de iqvoo 

' : ' : vyTluete@ 

! if f ticle i jane 3 oye 





3.6 


e Ri- 0.0 


a6) 0 -0.05<Ri<-0.15 


a Ri--0.3 


4 Ri:-0.4 
« Ri:-1.0 


* 4eoe 


Lne 





Figure 30. Effect of Stability Correction for 0637-0652 
21 Sept 73. 
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Figure 31. Non-dimensional Dissipation Rate as a 
Function of z/L (after Garratt). 
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(1974) described the influence of wind-wave coupling as 
having the same effect as would be due to the influence 
of stable stratification. The wave-related effect de- 
scribed by Davidson was a reduction in momentum transfer 
from that associated with neutral conditions. The re- 
sults of the present study are consistent with Davidson's 
findings in that there is a decreased momentum transfer, 


u as the surface is approached. If this were due to the 


«x! 
influence of wind-wave coupling, it would yield the ob- 


served deviation from the neutral slope, shown in Figures 


22 and 23. 


F. VARIATION OF TEMPERATURE FLUCTUATIONS 

Figure 32 represents the only temperature spectrum 
obtained from the five data periods which displayed the 
predicted -5/3 slope. Figure 33 is an example of the 
temperature spectra obtained from the majority of data 
runs. The paucity of temperature data precluded serious 


evaluations of the temperature fluctuation expressions. 
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Figure 32. Ten Minute Temperature Spectrum for 2323-2341 
20 Sept 73. 
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Figure 33, Ten Minute Temperature Spectrum for 2008-2018 
20 Sept 73. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The results obtained from the one poimitt spectral 
method are generally lower than those obtained by the 
intercept method. However, the simplicity of the one 
point spectral method would permit amore rapid deter- 
Mination of the dissipation and if a general estimate 
were ail that was desired, this method would be accept- 
able. The results of the intercept method show excel- 
lent agreement among levels and with results from other 
studies. The longer time averages (~20 min.) appear 
to be necessary to avoid acceleration effects at the 
lower level. 

The close agreement between dissipation values de- 
rived from the spectral estimates and the inner scale 
calculations confirm the spectral estimates and also the 
choice of the constant, a = 0.48. Results on the varia- 
tion of € with height also verified the formulations 
used in calculating the dissipation and flux. In addi- 
tion, the sensitivity of the results to the varying con- 
ditions of hydrostatic stability and, perhaps, wind-wave 
coupling support the predictions proposed by Deacon and 


Davidson. 
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B. RECOMMENDATIONS FOR FUTURE EXPERIMENTS 

It is strongly recommended that a real time spectrum 
analyzer be used in future observational experiments. In 
addition to the capability for instantaneous evaluation 
of spectral characteristics, the resulting spectral plots 
could be used for general estimates of the dissipation. 
If observational experiments are to be made on a more fre- 
quent basis, there is insufficient time to properly 
evaluate the data before the next experiment. The use of 
real time analysis techniques would preclude the need for 
digitization and tape manipulations and save valuable and 


costly computer time. 
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APPENDIX B 


A. CALIBRATION AND SCALING PROCEDURES 
1. Velocity Calibration 

Calibration of the velocity sensors was accom- 
plished using a manometer and a TSI Model 1125 Calibrator. 
The calibration procedure was to record woltage output ver- 
sus manometer height over a certain range of velocities. 
The computer program included in this Appendix was eiiees 
to generate the calibration curve for further use in estab- 
lishing the calibration factor. This procedure was followed 
for each of the velocity systems and resulted in the cali- 
bration curves included in this Appendix. 

The calibration curve obtained from this program 
resulted ina plot of v2 versus Yu where V is voltage and 


u is wind speed. 






= b 


wee Intercept 
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This curve describes a line 


2 


Vv = Ata) oCS 


+ b (1) 


For use in data analysis it was necessary to obtain a 


scale such that 


G& x ie" =) ut 
where C = calibration factor; e'.= voltage fluctuation 
and uu’ = velocity fluctuation. Differentiation of Equa- 


tion (1) yields 


2v dv = 1/2 eee Gia 

or 
2 
pew Vu edu au 
a 

or 

avul/? 

[ ] e' = y! 

a 


The calibration factor is then 


avul/? 


Cc a Umea ppm 


where "C" scales between volts and cm/sec. 
The procedures for actual calculation of this 


factor consisted of (1) obtaining the slope of the cali- 


bration curve (2) selecting a "u" valve equal to that 
observed during the measurement period and (3) entering 
the curve with this “u" value to obtain the corresponding 


voltage. 
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2. Scaling Procedures for Estimation of Epsilon 


In terms of the output of the spectral analyses 
programs the units associated with f$¢(f) or ko(k) were 
Kees). An examination of Equation (1), Section II, 
shows that the units should be am? SOS - Recalling that 
the calibration factor previously discussed scales between 
volts and cm/sec, it is necessary to multiply the € value 
by 


2 2 3/2 
(cm/sec ) _ |em/sec 


(memes) = 


to obtain the value of € in Eni seer L 

To obtain the correct value of € it is also neces- 
sary to consider the gains used during the record and 
analog-to-digital conversion phases. In addition a scale 
factor is needed that scales between the analog voltage 
and the digital output. A 10 db gain was used when the 
raw signals Were recorded and an additional 50 db gain was 
employed during the A-D process. The scaling between the 
voltage input and the corresponding digital output is ob- 
tained by considering that the maximum voltage input to 
the COMCOR Ci 5000 is 100 volts and this corresponds to 
a digital maximum of Dee on the XDX 9300. 

The total gain and scaling for the observational 


data then becomes 


digital a Value obtained 
WiWOLeEs) 2 Che G2 Fel volts ) from IBM 360 
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where V= 


To obtain the 


Signal from velocity unit 

amplification prior to record (10 db) 
amplification during A-D conversion (50 db) 
analog-to-digital scaling factor (273100) 


value which corresponds to that which can 


be used with the calibration factor previously obtained 


Vie = 


(IBM 360 value x 1/SC x 1/Gl x 1/G2) 


The final scaling then becomes 


3 
[iscare x Gainl x Gain2] x c| 


or in this instance 


23 3 
fs x O.1 x 0.02] x c | 
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Figure El. Real Time Spectral Analysis of Velocity Signal 
for 2008-2018 20 Sept 73. 
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Figure E2. 
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Real Time Spectral Analysis of Velocity Signal 
for 2023-2057 20 Sept 73. 
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APPENDIX F 


A. TAPEOUT PROGRAM 

TAPEOUT is a general purpose program used to determine 
the label and file contents of magnetic tapes. The program 
can take an otherwise unreadable tape and provide the user 
with Baer acd one information to permit access to the tape. 
TAPEOUT Ginn print a file number, density, parity, length 
in bytes of the record and a dump of the record. In this 
study the two key outputs were the record length and record 
dump. These outputs permitted the user to obtain the block- 
size and ascertain the location of erroneous records. The 
following is an example of the TAPEOUT deck set-up used in 
this study. 

//TJOHNSTON JOB (2413,0392,XM24), ‘JOHNSTON SMC 2561' 

//EXEC TAPEOUT, PARM='2.0,,2' 

//TAPEIN DD UNIT=2400-1,LABEL=(,BLP), VOL=SER=KARLYN 

/* 

TAPEOUT is part of the IBM 360/67 subroutine library 
and a discussion of the parameters (PARM=) is available in 
the computer center consultant's office and Anes in Com- 


puter Center Technical Note No. 0211-08. 


B. CONVERT PROGRAM 
The seven-track octal tape is converted into nine-track 
hexidecimal data samples usin¢e the CONVERT program. The 


block characteristics of the data samples are still 
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maintained on the nine-track tape but two additional para- 
Meters are affixed to each block. These two parameters, 
NCHAN and KMAX, are affixed to the beginming of each block. 
NCHAN is the number of channels of analog data digitized 
on the seven-track tape. KMAX is a numerical value set 
equal to the maximum number of samples per block. The in- 
sertion of these numerical parameters at the beginning of 
each block and the change in the block format by the FORM 
subroutine will put the seven-track tape in the proper for- 
mat for use in the UBCFTOR program. The CONVERT program 
used in this study along with sample JCL cards is given 


below. 


//JOHNSTON JOB (2413,0392,XM24),'JOHNSTON SMC 2561' 
//CONVERT EXEC FORTCLG, REGION. GO=100K 
//FORT.SYSIN DD 
DIMENSION IDAT(1024), DAT(1024) 
FACTOR = l. 
REWIND 2 
REWIND 4 
NCHAN = 8 
NRECL = 1024 
J = 0 
10 READ (2,20,END=60,ERR=40) IDAT 
20 FORMAT (20(50A4) ,24A4) 
J=J+1 
CALL FORM(IDAT,NRECL) 
DO 30 I=1,NRECL 
WRITE (4) NRECL,NCHAN,DAT 
GO TO 10 
40 WRITE (6,50) J 
50 FORMAT ('O',5X,'READ ERROR,RECORD NO.=',14) 
COM TO! 10 
60 WRITE (6,70) J 
70 FORMAT ('0O',5X,'END OF TAPE,RECORD NO.=',I4) 
END FILE 4 
RETURN 
END 


J* 


105 


“86% isneitibbs awit aid oqe2. Avkxe-entn of? wo Baaiaing 


, 2 7O4S30E THO w? 4 ey ~-#aGid. ws } + Bexitts site asad 


i 


-Mpalid Mose. to, painatpsd seit ot Bsxitin o2n ,RARR BoE 


Patidcipib 628o pollens Strats soungvea e727 ef 


z 
o 
dnd 
7 
o 
& 
i 


Mat peta THSEVVO! 4 i3'] f : ec) { ni ou 293. 





//GO.FTO2F001 DD UNIT=2400-1, VOL=SER=KARLYN,LABEL=(,NL), 

// DISP=OLD, DCB= (DEN=2, RECFM=VS, BLKSIZE=8192) 

//GO.FTO4FO01 DD UNIT=2400, VOL=SER=NPS170, DSNAME=DEREK, 

LABEL=(1,SL) ,DISP=(NEW,KEEP) , DCB= (DEN=2,RECFM=VS, 

// BLKSIZE=8200) 

It is noted that the user inputs to the program consist 
of the dimension of the arrays IDAT and DAT and the value 
of NRECL. IDAT and DAT should be dimensioned at least NSAMP 
x NCHAN thus signifying the need for recording the digitiz- 
ing parameters. The argument NRECL must be assigned a 
value equal to NSAMP x NCHAN. The value of the scaling fac- 
tor has been set equal to one in this study since it was 
found to be easier to perform all calibrations and scalings 
in one step following the use of the spectral analysis pro- 


grams. Computer Center Technical Note No. 0211-08 provides 


additional information concerning the CONVERT program. 


C. DIGITAL TAPE OPERATIONS 

Unlike the Hybrid computer system, where all tape mani- 
pulations were performed by the user with technical assist- 
ance from the computer enc uricy staff, the IBM 360 computer 
system manipulations were performed by the staff members. 
In order to have the tapes handled as desired, explicit in- 
structions had to be given. These instructions were given 
in two ways. The first instruction consisted of a Job Re- 
quest Form addressed to the computer center staff identify- 
ing the tapes to be utilized and whether they were to be 


read off of (RINGOUT) or written on (RINGIN). The second 
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instructional procedure was conducted usimg job control 
language (JCL) addressed to the computer and read as part 
of the submitted program. 

1. Job Control Language 

Job control language cards provide information to 
the computer about tape identification, tape disposition, 
tape data files, tape mount neaieres and data format. The 
cards are appended to the program and read through the "hot 
card reader" as part of the program, whereas the Job Request 
Form is submitted "over the counter" to the computer center 
personnel. An example of a JCL card for ttape processing 
would be: 

//GO.FTO8FOO1 DD UNIT=2400-1,VOL=SER=KARLYN, LABEL=(1,SL), 

// DSNAME=DEREK,DISP=(NEW,KEEP), DCB=(DEN=2,RECFM=VS, 

// BLKSIZE=8192) 

While the majority of the JCL eared remaims unchanged, there 
are portions which are subject to change from tape to tape. 
These portions are discussed below. 

FTO8 - This group indicates the tape mount on which the 
tape is to be placed. The two digit number (08 in this case) 
must match the number provided in the pragram's "READ" or 
"WRITE" statement. 

FOO1 - This group indicates the sequential number of 
puieees through the tape which has occurred up to this point. 
In this case it is the first pass. 

UNIT=2400-1 - This group specifies the tape as a seven- 


track tape. The nine-track designation is 2400: 
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VOL=SER=KARLYN - This group identifies the tape to 
be processed. In this case the name "KARLYN" signifies 
the tape is an "external" tape with the mame "KARLYN" 
physically appearing on the tape. An "imternal" tape 
would have a serial number such as NPS211. These tapes 
are kept on file in the computer center and cannot be 
Physically removed from the facility. 

LABEL=(1,SL) - This group identifies the tape file to 
be processed and specifies that the tape has a standard 
label. "External" tapes may be processed as no label 
tapes (NL) but "internal" tapes will be given a standard 
label before the user gains access to them. This is accom- 
plished by the computer center staff processing a "canned" 
program which places a standard label om the serialized 
tape. 

DSNAME=DEREK - This group specifies the data set name 
(DSNAME) applicable to the specified LABEL group. The 
data set name may be a unique name limited to six alpha- 
numeric characters. Once a file has beem given a data set 
name access to the data file can only be obtained if this 
same DSNAME is used. This process can be circumvented by 
the use of by-pass label processing (BLP) but does provide 
protection against accidental erasures. 

DISP=(NEW,KEEP) - This group specifies the tape dis- 
position once it has been processed. 

DCB= (DEN=2, RECFM=VS, BLKSIZE=8192) - This group is re- 


ferred to as the data control block. The density (DEN=2), 
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tape record format (RECFM=VS) and block size (BLKSIZE= 
8192) are stipulated using this group. The density speci- 
fies the number of bytes per inch, in this case the 

number 2 denotes 556 bytes per inch. This corresponds to 
a manual setting of the DENSITY dial on the face of the 
seven-track tape unit. RECFM=VS specifies that the record 
format is variable spaced, the result of recording and 
digitizing on the seven track unit. The block size speci- 
fies the number of bytes per block of data. If 2048 words 


per block were used, four bytes per word would result in 


8192 bytes per block. 
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APPENDIX G 


A. GENERAL 

All the UBC time series analysis programs (FTOR, SCOR, 
FC PLOT) were originally listed on faculty tape NPS216, 
file 1-3 respectively. fThe pagans have since been placed 
on computer center disk, MARY, and are stored in machine 
language which disallows a LAgehag of the program for user 
study. However, source decks for UBCFTOR and UBCSCOR have 
been retained by Dr. Kenneth Davidson of the Meteorology 
Department and by Dr. Noel Boston of the Oceanography De- 
partment. These source decks can be used to obtain a 
listing of the programs for individual study. The most 
recent computer center documentation on the use of the time 
series analysis package is given in the following memoran- 
dum written by Miss Sharon Raney of the computer center 
staff. This memorandum provides the information and JCL 
necessary for utilization of the disk-stored analysis 


programs. 


B. FTOR 

The preamble to the FTOR program is given following 
this discussion. It provides the user with the informa- 
tion needed to prepare the input data cards necessary for 
program execution. One of the most misleading parameters 


called for in the input values is the variable NBLOCK. 
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This value is provided by the user on the second data 
card with the value supposedly entered im columns 4-5, 
In actuality the value is read as a five digit integer 
thereby making the maximum value of NBLOCK 99999 vice 
99 as implied by the program preface. 

An important consideration involves the value of the 
variable MPRINT. This variable controls the computer 
printout of the FTOR program. If the user desires to sup- 
press the printout the insertion of a specific value can 
save considerable paper and cut computer time see or 
In the unsuppressed mode, the output will follow the format 
shown in Figure 13, Section IV. The suppressed output will 
consist of one line outputs for each data block which will 
inform the user that "coefficients for block __ have been 
Wiainttenw One en riko mu Ot the Out pUc. Capes. By uSiing ie hie 
suppressed output, computer turn-around time was signifi- 


cantly reduced. 


Cc. SCOR 

The preamble to the SCOR program also provides the user 
with the information necessary to compile the data input 
cards. The SCOR program requires the input of several 
variables whose values are determined during the A-D con- 
version process. One of the values on the first data card 
is the number of data blocks on which the spectral analysis 
will be performed. This value is set by the user and can- 


not exceed the number of data blocks processed by the FTOR 
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program. As was the case in the FTOR program the maximum 
number of blocks which can be handled is 99999 and not 99 
as the preface would lead the user to believe. 

As was discussed earlier in this thesis the output of 
the SCOR program consists of a computer printout anda 
graphic plot. The graphic plot consists of an output in 
which spectral density is plotted against frequency for 
all channels being processed. Since the original formula- 
tion of the SCOR progran, the computer facility at the 
Naval PAL eeetanat School has installed a new plotting 
routine called NEWPLOT. This routine limits the user to 
30 inches of plotting paper. The author found it necessary 
to submit a job request asking that plotting routine CAL- 
PLOT be used to prevent early termination of spectral 
plots consisting of four or more channels. The computer 
center staff has been alerted to this problem and any 
future user should consult with a staff member before 


using the SCOR program. 
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8 June 1973 


MEMORANDUM 


From: Sharon D. Raney 
To: Users of UBC Time Series Analysis Package 


Subj: New Operating Instructions for UBC Time Series Analysis Package 
1. The new FTOR deck setup is given in Figure l. 


The dimension size of array DATA must correspond to the block size 

of the data on the input tape. This tape is generally the data as 

converted from the seven track XDS generated tape to the nine track 
IBM/360 tape. 


Formerly, this array was dimensioned in subroutine OCEAN2 and had 
to be changed whenever the sample size changed. This meant the 
entire source deck for FTOR had to be used for each run. 


FTOR is now compiled in FORTRAN IV (H) which generates a more 
efficient code and consequently speeds up execution time. 


2. When the same dimension size is to be used repeatedly, Figure 2 
shows how to obtain an object deck and how to execute FTOR with the 
object deck. 


3. Figure 3 shows the deck setup for SCOR. 


4. Notice the change in the format of the Job Control Language (JCL) 
to define the tape files. These changes enable more efficient 
processing by the operator. Note the use of the IN parameter in 
the LABEL parameter. This must be specified when the file is to 
be read only. OUT is used when the file is to be written only. 
All data sets are PASSed not KEEP. The UNIT-AFF is used on sub- 
sequent files of a tape as is the VOL=REF. 


5. Any questions concerning these new operating instructions should 
be directed to Sharon Raney at X2631 or Ingersoll Hall Room 107. 
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Figure G2. 
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FACTOR PROVIDES THE DIMENSIONAL SCALING 
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momentum transfer during periods of stable stratifications were consistent 
with wind-wave coupling effects described by Davidson. The shape of a 
spectrum of temperature fluctuations was in agreement with predictions. 
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